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AN EXPERIMENTAL AND COMPUTATIONAL INVESTICATION OF THE R.OW FIEU) ABOUT A TRANSONIC AIRFOIL 
IN SUPERCRITICAL FLOW WITH TURBULENT BOUNDAR"-LAYER SEPARATION 


Morris W. Rubesln,* Arthur F, Okuno,** Lionel L. Lew, .Ir.,** 
John B. McDevitt,** and H. Lee Seegnlller** 

Ames Research Center, NASA, Moffett Field, California S4035 


Abstract 

A combined experimental and computational 
research program for testing and guiding turbulence 
modeling within regions of separation Induced by 
shock waves Incident on turbulent boundary layers is 
described. Specifically, studies are made of the 
separated flow over the rear portion of an 18X-thlck 
circular-arc airfoil at zero angle of attack in high 
Reynolds number supercritical flow. The measure- 
ments include distributions of surface static pres- 
sure and local skin friction. The Instruments 
emnloyed Include high-frequency response pressure 
cells and a large array of surface hot-wire skin- 
friction gages. Computations at the experimental 
flow conditions are made using time-dependent solu- 
tions of ensemble-averaged Navler-Stokes equations, 
plus additional equations for the turbulence 
modeling. 


X • distance along chord from leading edge 
Ax “ effective length of skin-friction gage 
y » distance from surface normal to chord 
“ location of separation streamline 
6 ” thickness of viscous region 

6* • incomnressible displacement thickness (Eq. (5)) 

i 

c « eddv diffuslvltv 

u “ viscosity at the airfoil surface 

w 

V « kinematic viscosity in the free stream 
« free-strean density 


Nomenclature 

A ■ van Driest damping length 

Cp " pressure coefficient, (p^ - P^)/q„ 

c “ airfoil chord length 

c, • local skin-friction coefficient based on the 
free-stream dynamic pressure 

I • wire-gage current 

k ■ roughness height 

t • mixing length of turbulence 

M “ free-str«»am Mach number 

Pr ■ Prandtl number 

p^ • local surface pressure 

p • free-stream pressure 

2 

q^ ■ free-stream dynamic pressure, (l/2)p^u_^ 

R “ resistance of wire gage 

Re “ Reynolds number, u c/v 

c,» •' • » » 

RSg <• Reynolds number based on boundarv-layer 
momentum '-hlckness 

R^ ■ resistance of wire gage at 0°C 

u • local velocity parallel to chord 

u, • velocity at edge of viscous region 
0 

u “ free-stream velocity 

Oi 

V “ local velocity normal to chord 

*Scr1or Staff Scientist. 

♦♦Research Scientist. 


p ■■ density at the airfoil surface 
w 

T • surface shear stress 

w 

T • surface shear stress on surface with constant 
wo 

pressure 

Introduction 

An understanding of the fluid-dynamic mech- 
anisms controlling the separation of boundary 
layers from aerodynamic surfaces is of ma1or Impor- 
tance In aeronautics. The forces and moments on 
aerodynamic bodies at angle of attaca and the per- 
formance of diffusers within wind tunnels or air- 
craft engines are often governed by the existence 
and location of separation. Thus, the oredlction 
of the position of the onset, the behavior of the 
bound.arv-layer separation, and the reattachment 
processes have occupied the fluid mechanics commun- 
ity for many years (e.g.. Ref. 11. Until recently, 
however, the techniques of predicting separation 
have been quite approximate, at best involving itera- 
tion between essentially Invlscid and viscous 
regions of the flow field. ^ In this Iteration pro- 
cess the viscous region is tightly coupled to the 
inviscld region, whereas the reverse Is not true. 

This unequal sensitivity had led to numerical diffi- 
culties, and the iterative process is not generally 
successful. 

Recent developments in computer technology and 
numerical analysis, however, have made it possible 
to avoid this iterative process by permitting the 
time-dependent solution of tlie full Navler-Stokes 
equations for the entire flow field. ^ For two- 
dimensional laminar boundary lavers, these solutions 
have been found to be quite accurate. '* When atten- 
tion is turned toward turbulent boundary lavers with 
eddy viscosity representations of the Reynolds 
stress terms in the ensemble-averaged Navler-Stokes 
equations, it Is found that these empirical methods, 
so useful for equilibrium boundary layers, do not 
work well when separation occurs. 


i 

I 


In an effort to improve turbulence modellnK for 
boundary layers scparatlnR at sonic flight speeds, 
the Ames Research Center has initiated concomitant 
proftrams in experimental and numerical fluid dynam- 
ics in which the experimental results are used to 
help guide the development of turbulence models. 

The first part of this program was reported in 
Ref. 7. The experiment involved the measurement of 
the static pressure distribution over the surface of 
an 18X-thick circular-arc airfoil in the free-stream 
Mach number range from 0.71 to 0.79 for a range of 
Reynolds numbers. These measurements were supple- 
mented by flow visualization irith surface oil 
streaks and shadowgraphs. The airfoil was tested in 
the Ames High Reynolds Number Channel, a facility 
that has top and bottom walls contoured to the shape 
of the streamlines computed to exist one-half the 
channel height away from the chordline of the air- 
foil in an invlscid flow for a specified design Mach 
number. The Reynolds numbers employed in the experi- 
ment were as high as 17 x lO^ , based on free-stream 
flow conditions and the chord length. At the high- 
est Mach number tested, the standing shock waves on 
the airfoil only extended two-thirds of the distance 
to the contoured walls so that wind-tunnel blockage 
was not a factor in these experlnents. In parallel 
with the experiments, computations of the flow field 
about the airfoil were performed through time- 
dependent solutions of ensemble-averaged compressible 
Navier-Stokes equations in finite-difference fom. 

The turbulence models used in these calculations 
Included the usual algebraic formulations for eddy 
viscosity,*® which imply the turbulence is in equilib- 
rium »»lth the mean flow, and semi-empirical turbu- 
lence lag models that attempt to account for turbu- 
lence history. In addition, some unique eddy viscos- 
ity formulations were required in the separated 
region downstream of strong standing shock waves; 
these are presented in the present paper in the sec- 
tion devoted to the flow-field computations. 

In Ref. 7 detailed comparisons were made between 
the measured and computed surface pressure distribu- 
tions for two Mach numbers at a chord Reynolds num- 
ber of 4 X 10®. In either case, a standing shock 
wave terminates the supercritical flow, that is, the 
supersonic flow region that develops adjacent to the 
airfoil. At M • 0.74, the standing shock wave is 
too weak to separate the bimndary layer, and separa- 
tion only occurs very close to the trailing edge 
where the pressure rises quite rapidly. The compu- 
tations, with the different eddy viscosity models of 
turbulence, predicted shock-wave locations about 
5-10* of the chord downstream of the actual location. 
Modifications of the turbulence model that improved 
the prediction of the shock-wave location resulted 
In too low a pressuie recoverv at the trailing edge. 

At M >0.79, the standing shock wave Induced sepa- 
ration at about x/c 0.675. When the turbulence 
models were adjusted in an attempt to fix the stand- 
ing shock wave where the experimental wave is 
located, the predicted pressure recovery in the 
separated region was much too high. 

Surface pressure distributions, while contribut- 
ing to an overall assessment of the accuracy of a 
specific turbulence model computation for a given 
flow condition, contain insufficient information to 
guide Improvements necessary In the turbulence model. 
Traditionally, this guidance has been provided bv 
detailed mean velocity survevs through the boundarw 
layer (e.g. , Refs. 10 through 14), with the nensiirc • 
ments close to the surface being particularly useful. 


Because the current experiment Involved a very snal 1 
model operating at a high unit Reynolds number 
(80 X 10® /n), the boundary layers on the model were 
too thin (0.4 nm) to survey with impact pressure 
probes, hot wires, or a laser doppler veloclmeter In 
either the viscous sublayer or the fullv turbulent 
wall region of the boundary laver. i^urther, intrus- 
ive probes, such as Impact pressure probes or hot 
wires, together with their support mechanisms, were 
considered possible sources of flow Interference at 
the transonic flow conditions of the present 
experiment. To avoid such interference and, yet, to 
gain insight into the mechanisms taking place within 
the boundary laver close to the surface, reliance 
was placed on measurements of the local skin 
friction. 

The flush heated wire gage described in Ref. 15 
was particularly suited to the task of measuring the 
local skin friction because of Its compact configu- 
ration and insensitivity to streanwlse pressure 
gradients that wojild occur on the airfoil. A design 
decision had to be made as to whether the individual 
gages would be (1) mounted on buttons that could be 
calibrated in channel flows and then placed in cavi- 
ties on the wing surface; or (2) mounted directly on 
the surface of the wing and calibrated indirectly by 
operating the wing under conditions where boundarv- 
laver calculations of local skin friction were 
expected to be accurate. The latter option was 
adopted as the simpler approach, and one which could 
yield an aerodynamically smoother surface. Thus, in 
effect, the gages read akin-friction values relative 
to those that existed on the wing at free-stream 
Mach numbers where the external and boundary-layer 
flows were believed to be understood. 

Since all the Instrumentation employed in the 
present investigation had response times of less 
than 0.01 sec, dynamic as well as static measure- 
ments were taken to determine if shock-induced sepa- 
rations that appear to be steady were really so. 
Thus, dynamic measurements of skin friction are 
presented. 

Because the experiment described herein uti- 
lized the same facility reported thoroughly in Ref. 
7, the experimental emphasis here will be on the 
skln-f rlrtion model. For completeness, the turbu- 
lence model used In the calculations presented in 
this report is given In detail because It forms the 
basis for some discussion of the results. Finally, 
the dynamic behavior of the surface measurements Is 
discussed. 

Experimental Procedures 
Ames High Reynolds Number Test Channel 

The Ames High Reynolds Number Test Channel was 
used for the skln-f rlctlon experiments described. 

The channel consists of a large high-pressure air 
supply, a test section, and a large-volume vacuum 
vessel, described in detail In Ref. 7. A schematic 
diagram of the test section of the channel Is shown 
In Fig. 1. The air enters the test section through 
a hellmouth, passes through the test section and a 
speed control region before entering the diffuser 
that leads to the vacuum vessels. For the current 
experiment, the pressure distribution and skln- 
frlctlon airfoil models spanned the test section and 
were mounted from the side walls. The upper and 
lower test section walls were contoured as described 
in the previous section. In addition to the 


REPRODUClBn.ITY OF THF 
ORMNAL PAGE 18 P(K)R 


r 



i 


J 1 L 

I 0 I 



Pig. 1 Schematic diagram of Ames High Reynolds 
Number Test Channel 



Fig. 2 Skin-friction model instrumentation 


contouring, the upper and lower walls of the test 
section could be pivoted about their leading edges to 
permit compensation for the displacement effects 
created by the boundary-layer growth on the four 
walls of the test section. The Mach number in the 
channel was controlled by adjusting the flow area in 
the speed control section with Inserts in the top 
and bottom walls and by translating a wedge axially 
along the test section centerline. 

The air and vacuum supplies of the Ames High 
Reynolds Number Channel are adequate to permit opera- 
tion at a relatively high unit Reynolds number of 
85 X 10*^/m for a maximum period of 45 sec. During 
this time interval, Joule-Thompson cooling of the 
air supply causes the stagnation temperature of the 
facility to drop by about 10®C. 

Experimental Models 

Two models having identical external configura- 
tions were utilized in the experiment. Each is a 
circular-arc airfoil, 18Z thick, with a chord length 
of 0.2012 m. The first model, designed for measure- 
ments of the pressure distribution about the model, 
was made of steel with a surface finish of 0.000813 
mm. The model contained 47 pressure taps on one sur- 
face and 3 on the other side to allow aligning the 
alrloll to an angle of attack of 0°. Details of 
this model and of the pressure distribution experi- 
ment are given in Ref. 7. The second model, used 
for this skin-friction experiment, contained 7 pres- 
sure taps in addition to an array of 30 hot-wire 
skin-friction gages set into blocks of plastic. The 
purpose of the pressure taps was to ensure an align- 
ment of the model to a 0*^ angle of attack and to 
verify that the pressure distribution on the skin- 
friction model, at a few discrete stations, agreed 
with the pressure measurements at the same stations 
obtained on the pressure distribution model for the 
same free-stream conditions and top and bottom test 
section wall settings. Although the models spanned 
the test section of the channel, the flow field in 
the vicinity of the center span of the model away 
from the side-wall boundary layers was found to be 
two-dimensional, even under conditions where large 
regions of separation occurred over the airfoil and 
on the side walls. ^ 

The schematic drawing of the skin-friction 
model, presented in Fig. 2, shows the sections of 
two rather large cavities in the model (shaded areas) 
containing polystyrene plastic (rexolite) that was 
molded into place and which supports the lead-in 
wires and the skin-friction gage wire elements. 


Table 1 Location of skin-friction gages 


Wire 

x/c 

Wire 

x/c 

Pressure 

x/c 

no. 

no. 

orifices 

1 

0.276 

16 

0.651 

1.2 

0.1 

2 

0.292 

17 

0.676 

3 

0.3 

3 

0.308 

18 

0.713 

4 

0.5 

4 

0.324 

19 

0.744 

5 

0.7 

5 

0.437 

20 

0.775 

6,7 

0.9 

6 

0.453 

21 

0.805 



7 

0.469 

22 

0.821 



8 

0.484 

23 

0.836 



9 

0.500 

24 

0.851 



10 

0.516 

25 

0.866 



11 

0.531 

2C 

0.881 



12 

0.547 

27 

0.896 



13 

0.562 

28 

0.911 



14 

0.589 

29 

0.926 



15 

0.621 

30 

0.941 




Photographs of the model are shown in Figs. 3 and 4. 
The dull spanwise streaks are filled cavities that 
contain tubes leading to the pressure taps. The 
wires composing the skin-friction gages were placed 
on the surface of the plastic Inserts in a direction 
normal to the stream between pairs of lead-in posts, 
and epoxy plastic was molded in between the wires as 
shown in the inset on Fig. 2. This technique for 
construction of the skin-friction gages is identical 
with that developed in Ref. 15. On the current 
model the gage wires were made of platinum having a 
diameter of 0.0254 mm and lengths in the spanwise 
direction of about 6.35 mm. One problem with this 
type of model construction is the maintenance of an 
aerodynamlcally smooth surface at the Joints between 
the polystyrene and steel, and in the epoxy between 
the individual wire gages. Because of the relative 
softness of the plastic surfaces and the rather 
delicate nature of the wire gagi-s, it was decided 
to make surface contour measurements without actu- 
ally touching the surface. A low-intensity, but 
narrow laser light beam was projected on the surface 
of the mi.'del and the position of the reflection was 
noted on a screen across the laboratorv. The local 
surface angle of the model was determined from the 
geometric arrangement of the laser, model, and 
screen. From these measurements it was a simple 
matter to deduce surface wavlness or roughness 
heights of about 0.0004 mm. It was found that the 
model had a waviness of less than 0.004 mm in an 
average distance of 3.18 mm between the wire gages 
(see inset in Fig. 2). The step heights in the 



Fig. 3 Photograph of skin friction codel. 


The buried-wlre skln-f rlctlon gages used In 
this experiment are still In the development stage; 
therefore, their calibration was conductsd as an 
intimate part of the main experiment. Details of 
the calibration procedure and related discussion are 
given in the results section. 

Computational Method 

The computations of the flow field about the 
airfoil at • 0.788, where boundary-layer separa- 
tion is Induced by the shock wave, were performed 
using a form of Che computer program described in 
Ref. 9. The computer program utilizes an explicit 
finite-difference method to solve the time-dependent, 
two-dimensional ensemble-averaged Navier-Stokes equa- 
tions applicable to compressible turbulent flows. 

The turbulence is modeled with an algebraically 
expressed eddy dlffuslvity model, termed Model 3 in 
Ref. 9. This model assumes that Che turbulence is 
in equilibrium with the local mean flow. With ref- 
erence to the following sketch, the eddy dlffuslvity 
expresrlons for the various regions indicated are as 

follows: 

( 



Fig. 4 Detailed photograph of skin friction model 
showing lead-in posts and gage wires. 


surfaces where the plastic and steel met were less 
than 0.0003 mm. For reference, the criterion for 
aerod ynam ic smoothness in an attached boundary layer 
is * c j/2 (u_k/v^) < 5, where k is the allowable 
surface roughness height.' In the regions of 
attached flow for the conditions of this experiment, 
the values of k are about 0.0023 mm or one-half 
the actual surface waviness. It can be expected, 
from the data summarized in Ref. 1, that the surface 
waviness on the skln-f rlcti'— .»odel could cause the 
measured skin friction to bt' about SZ high in the 
regi PS of attached flow. For separated regions, 
the surface waviness should be less influential. 
These are acceptable errors for an heuristic experi- 
ment such as this. 

For future reference, the numbering and corre- 
sponding positions of the skin-friction gages are 
shown graphically in Fig. 2 and are also given in 
Table 1. 



Inner region I 


where 


r/_ \2 /. 


(1) 

0.41y(l - exp(-y/A)) 

(2) 

u.. 

, - 26 yp / T 

0 W W' 

(3) 


This is essentially the classical van Driest model 
for the inner portion of a turbulent boundary layer. 
The term 3v/3x has been added to Eq. (1) to 
account for departures from true boundary-layer flow 
as the point of separation is approached. The van 
Driest damping length A is expressed identically 
as on a flat plate, with no account being made for 
the effect of surface pressure gradients as is sug- 
gested, for example, in Ref. 10. 


II Outer region of boundary layer and wake 
0.0168U.6* 

£ m 6 1 

1 - l(y - yt,s)/«l« 


(4) 


where vq^ Is the location of the streamline bound- 
ary of the separated region and 4 is the overall 
thickness of the viscous region. The displacement 
thickness 4* is defined as 
1 


4* 

1 



(1 - u/Uj)dy 


(5) 





Equation (4) Is based on the Clauser model, but 
accounts for Intermlttency of the turbulence In the 
wake with the term contalnlnq y. The usual defini- 
tion of 4* Is modified with the lower limit of 
Integration moved from the surface to the separation 
streamline y^g. 

III Separation bubble wall region 

c - 0.0168u^4*{ (y/Vpg)[l - exp(-v/A)])^ (6) 

IV Separation bubble wake region 

E - 0.0168u^«* (7) 

Discussion of Results 
Surface Pressure Distribution 

The transonic flow over the circular-arc, 18Z- 
thlck airfoil was studied In Ref. 7 with the pres- 
sure distribution model for a rather extensive range 
of Mach and Reynolds numbers. In these experiments, 
the upper and lower walls of the channel were con- 
toured to conform to the free streamlines that would 
exist over this airfoil for free flight In an invls- 
cld fluid at ■ 0.775. The channel walls and 
their positioning «rere kept constant for all the 
Mach and Reynolds numbers of the tests. Mach number 
was controlled by adjusting the dimensions of a down- 
stream nozzle throat and the position of the trans- 
lating wedge (see Fig. 1). The Re^ntolds number was 
controlled by setting the stagnation pressure. The 
stagnation temperature Is uncontrolled, but reaialned 
close to ambient and was continuously monitored. 
Supercritical flow, where the local velocity at the 
surface of the airfoil exceeds the local sonic veloc- 
ity, was found to occur above - 0.71. Between 
” 0.76 and 0.78, an unsteady, asymmetric periodic 
flow occurred. At Mach numbers below the unsteady 
flow regime the boundary layer remained attached 
until the trailing edge was approached; in the 
unsteady regime, separation alternately switched 
from the trailing edge to the base of the shock wave; 
and above M > 0,78 separation was fixed at the base 
of the shock wave. 

Since the objective of the present experiment 
was to define the local skin friction for the case 
of nominally steady shock-induced separation. It 
was necessary to test at a Mach number greater than 
0.78. Thus, the skln-frlctlon experiments were con- 
ducted at M ■ 0.783 for chord Reynolds numbers 
ranging from 8 x 10^ to 14.3 « 10^. For calibration 
of the skln-frlctlon gages, the skln-frlctlon model 
was also run at M_ - 0.682 where the flow Is sub- 
critical and boundary-layer theory can be expected 
to apply over almost all of the airfoil with reason- 
able accuracy. Because the Mach numbers employed In 
these tests differ from that used to design the con- 
tours of the upper and lower walla, and also because 
viscous effects including the rather large separation 
zones would also cause the actual streamlines away 
from the airfoil to deviate from those based on 
Invlscld theory. It cannot be expected that the pres- 
sure and skln-frlctlon distributions measured In 
this experiment would conform closely to free-fllght 
conditions. It would be expected, however, that the 
behavior In the two types of flows, that Is, In the 
specific channel configuration and In free flight, 
would be sufficiently similar so that the separation 
mechanisms studied In one flow field would be repre- 
sentative of the other. In the channel, strict 


application of the computational method for calculat- 
ing the turbulence separation on the airfoil 
requires introduction of the actual upper and lower 
wall contours as boundary conditions In the computer 
program. This Is being done at present; however, 
the computer results given In this paper were based 
on free-fllght conditions and when cosq>arlng the 
data and the computations this Innate difference 
must be considered before drawing conclusions. 

Representstive surface-pressure distributions 
on the airfoil for the two test Mach numbers of this 
experiment are shown on Fig, S. The ordinate Is the 
pressure coefficient Cp and the abscissa Is the 
station represented as a fraction of the chord. The 
data for M " 0,784 and Re. " 10.7 « 10*’ Indicate 
shock-induced separation occurred at about 
x/c “ 0.675. Bevond this station the pressure 
recovery Is quite poor. The pressure measurements 
on the skln-frlctlon model (filled symbols) show 
that the flow conditions over the two models could 
be reproduced quite well. The subcrltlcal case is 
represented by the data at M • 0.682 and 
Re^ ^ ■ 12,4 * 10®. At this Mach number, signifi- 
cant pressure recovery occurs toward the trailing 
edge. Again, there is good agreement between the 
pressures measured on the oressure distribution and 
skln-frlctlon models (open and filled symbols, 
respectively). The data obtained at other Reynolds 
numbers (3,1 « 10® < Re^ ^ i 14.3 x 10®) are quite 
similar to those shown In Fig. 5. The positions of 
separation indicated on the figure for the two test 
cases were determined from oll-fllm photographs.^ 
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Fig. 5 Pressure distributions on airfoil for 
tralllng-edge and shock-induced 
separation 

Calibration of Skln-Frlctlon Pages 

The principle of the heated wire skln-frlctlon 
gage requires that Its dimension In the streamwlse 
direction he saull compared to the boundary-layer 
thickness. As the effective dimension of the gage 
Is composed of the actual dlaswter of the heated 
wire plus the region of the substrate surface whose 
temperature Is elevated because of the heat con- 
ducted upstream from the wire. It Is necessary for 
the gage wire diameter to be very fine and for the 
substrate to have a very low thermal conductivity 
(Refs. 15 and 16). IJhen these conditions are met, 
the heat transferred from the heated wire element 
forma a thermal boundary layer over the gage that 
lies within the viscous sublayer Immediately adja- 
cent to the surface. Because the flow here is 
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eusentiallv laminar, the heat-transfer nvchanihms 
arc governed by the fluid properticB and the veloe- 
lt> gradient normal to the surface, or skin friction, 
that exists iiist upstream of the gape. Tlius, the 
power lost from the gage, per unit of temperature 
rise over the unheated airfoil surface at the same 
station, becomes a measure of the local skin friction. 


The measuremt>nis of skin friction arc sensitive 
to the local surface-pressure gradients as wis sho\m 
In Refs. 15 and 17. Equation (5) of Ref, 15 con- 
tained a typographical error of one order of magni- 
tude foi the pressure gradient tern. To correct that 
error, and to put the equation for pressure gradient 
effects Into a form useful for the present experi- 
ment, the effect of pressure gradient is expressed 
as _ 

, 1''^ . '1C 

(^x/c) 


- 1 - 0.312 


(c,/2)Re^ Pr, 
. f c,~ J 


c, dfx/c) 


(«) 


where t„ Is the actual shear and Is the 

shear induced from the gage when the e?fects of sur- 
face pressure gradients are neglected, Iiicidentallv, 
when the equations of Ref. 17 were put In the form 
of q. (8), the constant in the picesure gradient 
term was 0.344. 


Equation (8) shows that a smaller 4x tends to 
reduce the 'ffect of nonzero pressure gradients, 
though the affect of Ax la only to the 1/3 power. 

It also sh- ws that at separation, whore Cr 0, 
pressure gradients dominate the value of the skin 
friction. To estiimate the effects of pressure grad- 
ients In the present experiment, Eq. (R) was applied 
to two portions of the airfoil. At M ■ 0.784, Eq. 
(81 was evaluated at x/c • 0.64, the location of 
the foot of the shock wave. With estimated values 
of cf “ 0,003 and Ax ■ 3 x (wire diameter),*^ It was 
found that the resulting was 2.2X lower than 
T^(,. the value of shear stress that would be inter- 
preted to exist when the effects of pressure grad- 
ient are Ignored. In the region after separation, 
the pressure gradient was much lower, but skin fric- 
tion magnitudes were also expected to be lower, for 
example, C£ 0.0003. Here errors Introduced by 
neglecting the effects of pressure gradient were 
about 0.7X. These estimated errors ^rero sufflclentlv 
small to lead to the decision to neglect surface- 
pressure gradient effects in the data reduction. 

The first step In calibrating the skln-f rlctlon 
gages was to determine their electrical resistance 
as a function of temperature. This is needed in 
order to use the gages as resistance elements as well 
as heater elements. Tills step was accomplished by 
placing the skln-frlctlon model In an environmental 
test chamber capable of malnlaininK temperatures 
between -130° and 300°C. Two thermocouples were 
located on the plastic Insert, one on the stec-1 air- 
foil, and one in the air Just above the plastic 
insert. A small fan Inside the chamber circulated 
the air to minimize temperature gradients. As a 
consequence, the thcrmccouples showed departures 
from tlielr mean temperature rending of only about 
0.06°C. The resistances of the wires were read and 
recorded cn the channel Instrumentation system 
(electronic mil llvol traeter plus magnetic tape 
recorders) so that ihe calibrations and application 
of the gages used the aane rend-nut pquipmepr. 

The results of the electrical resistance cali- 
bration of the 30 skln-frlctlon gages are shown on 
Elg. 6. For comparison, the resistance of n free 



Fig. 6 Temperature dependence of the electrical 
resistance of the skln-frlctlon gages 

wire loop of the sairn* platinum wire uaed for the 
gages is also shown. The referent e temperature wrs 
chosen to be 0°C, and measurements were made up to 
75°C. It Is observed tli.-'t the wire gages mounted on 
the airfoil model showed greater changes of resis- 
tance with temperature than did the free wire. This 
Is explained by the stretching of the mounted wires 
resulting fmm the oolvstvrene having a larger 
linear thermal expansion coefficient (7 ■ lfl"5/oo) 
than does the platinum (O.S** » 10~^/®C). Although 
the wires gave Individual readings of resistance 
change ulthln about IT ef each othor, the individual 
wire c.il Ihrafions were used in the data reduction to 
account for the small differences in resistance that 
existed between the wires, and which could cause 
tcnneratiire differences of a few degrees If n single 
mean curve were used. 

The skln-frlctlon model was then mounted in the 
Ames High Reynolds Humber Channel and the gage wires 
were electrically connected to be used either as 
resistance elements, for surface temperature mea- 
surements, or to be driven bv standard constant tem- 
perature mode hot-wl-e equipment, where the Instan- 
taneous voltage required to rainteln a preasslgned 
overheat is measured and recorded. As six DISA sys- 
tems were available, onlv six wires were heated 
during anv given run. This was not a serious limi- 
tation because It was shown In Ref. 15 that heating 
a gag.' immedlatelv upstream of another sage could 
aftect the latter's readings. For the measurement 
of the magnitude of skin friction, ■everal wires 
were left unheated between the heated ski n-f rlctlon 
gages to minimize the Interfcrt-iice. This inter- 
ference effect. In fact. Is the basis of the sepa- 
ration detecto' explained In Ref. 15 ind was used In 
the pn-sent experiment to define the sign of the 
skin-f r let Ion measurements In the separated zones. 
Unheated wires downwind of s heated gago indicated 
temperatures in excess of what would have been 
expected at that location In the ahaence of heating. 

As stated earlier, boundarv-laver cnmputatlot.s, 
under (.ondlt’ons where boundarv-laver theory is 
expected to give good Indications of the skin fric- 
tion, were used to "callhratc" the skln-frlctlon 
gages. To gain some Insight into the accuracy of 
this process, two different boundarv-laver programs 
containing different turbulence models were com- 
pared for the subcritlcal pressure distribution 
(M^ « 0.6fl2> ihown on Fig. 5. 

One pr-’g'-e-, H^tryin-Sheaffer, is an 

extension of Ref. 18 that contains a two-laver alge- 
braic eddv vlsconltv model similar to that of Ref. 

10 thereby, assumes that the turbulence Is 



always In equilibrium with the mt>an flow. The second 
computer program allows for departures from equilib- 
rium between the turbulence and the mean flow b’’ 
IntrodiielnR two rtlfierentlal equf.tlons for the Inten- 
sltv and scale oi rhe turbulence.*** This latter nro- 
gram Is an extension of the work of Saffman*® .nnd, 
for reference, is called the Modified Saffman model. 
In Its present form, it can alsr account for the 
streamline curvature along the surface of the air- 
foil as an option. 

The distributions of skin friction that resulted 
from applying these programs to the airfoil with a 
pressure distribution corresponding to « 0.ft82, 
Fig. 5, are shown In Fig. 7 for the case where the 
ch.innel was operated at a total picssure of 60 psin. 
The Harvln-Sheaf for program required the positions 
of the beginning and the end of transition to turbu- 
lence to be specified. The transition locations were 



Fig. 7 Calculated skin friction on airfoil for 

conditions of an attached boundary layer, 

M - 0.682, Re - 12.4 >• lO*’' 

- c,- 

estlmated from plots of the power lost per unit tem- 
perature rise, 1*R/4T, recorded by the skln-frlction 
gages. The location where the rR/A.T began rising 
rapidly with distance along the chord was chosen as 
the point of onset of transition. The location of 
tlie end of transition was taken where the I*R/4T 
reached a first maximum. Ifhen the skin-friction gage 
data at the four different stagnation pressures or 
Reynolds numbers were examined In this way. It was 
noted that transition onset could be represented 
reasonably well by a value of Re^ ■ 523 and the end 
of transition by Rc0 ■ 1567. Of course, not all the 
data agreed exactly with these criteria but the d'.ta 
could be bounded, for example, on the high side hv 
Re0 ■ 666 and Re^ •• 1699 for the onset and end of 
transition, respectively. The results from the 
Marvln-Sheaf fer program base«' on both acts of transi- 
tion criteria indicate In Fig. 7 that skin friction 
downstream of transition Is relatively InHensltlve 
to small changes In transition location. At the 
trailing edge the two calculations differ hv ahnut 
101, but It la not clear why differences caused bv 
the location of tranaitlon near the leading edge cf 
the airfoil wou’d be Inreeat fartheat downatrean 
from transition. 


To ahow how sensitive the predicted skin fric- 
tion Is to the narticulnr raoilel of turbulence 
emoloved in the cjlcul at ions , computations based on 
the Modified Snffman model sre nlsn .ihown on '^Ig. 7, 
Because the earlier Mar vln-Sheaf f er compuiallous did 
not ahow murh 8ensltl''ltv to the tranaitlon loca- 
tions, no particular care w.as t.-ikm to match the 
transition to the data. Tr.'frsltlon wss arbitrarily 
Initiated at x/c • 0.02 and allowed to run its 
courae within the computer program. One Modified 
Snffman curve Includes the effect of the convex sur- 
fnee curvature of the airfoil whereas the ot*’er 
treats the airfoil boundary laver as one on n flat 
plate experiencing the pressure distribution occur- 
ring on the airfoil. A comparison of these two 
curves Indicates that the effect of the convex sur- 
face curvature Is to reduce the local skin friction 
bv an amount that Increases with distance along the 
airfoil surface, or more directly, the boundary- 
laver thickness. At the trailing edge, where the 
boundary laver is the thickest, the difference 
caused bv the curvature la about 30t of the rela- 
tively small values of sheat occurring there. At 
x/c “ 0.7, the curvature effect Is onlv about 5X. 

Prom the curves shown in Fig. 7, It is 
observed that bevond x/c • 0,3, the region of 
interest in the senarated flow test, the results of 
the Marvln-Sheaf fer program and the Modified Saffman 
program with surface curvature agree to within lOZ. 
Thus, either program could be used for "calibrating" 
the aV.ln-f rlctlon gages; because the Varvln-Sheaf ftr 
program runs verv quicklv. It was used tor evaluat- 
ing the shear on the airfoil for channel total pres- 
sures of 15, 30, and 45 psla. In addition to the 60 
pala values shown In the figure. 

The skin-friction gage calibration ororesa con- 
sisted of testing the airfoil at M^ ■ 0.682 at a 
series of Reynolds numbers. During each test, aix 
widely separated skin-friction gages were heated and 
the power loss per unit of temperature rise (I^R/AT) 
within each heated gage was determined. The temper- 
ature of the heated gages and of those unhested 
wires measuring the airfoil surface temperature lust 
upstream of the heated gages (i.e., the AT) were 
determined from the electrical resistance R mea- 
sured for each wire and from the turves of Fig. 6. 




Fig, 8 Callbrat*oii mr-es ot representative skln- 
'r' flon gages 


For each successive run, another aeries of nix (taxes 
was heated, and so on until all the xaRes were heated 
I or a given Reynolds number. Five rupa w.-re needed 
nt each Reynolds number to cover the 10 I'.ages on the 
model and a sixth run was made to meanure llie air- 
foil Burfare temperature without heatintt anv of the 
wires. This process was repeated fur the four 
Reynolds numbers. When the measured power loss ner 
unit of temperature rise for each Rage was plotted 
against the calculated values of ^ for the 

corresponding test conditions and gage locations, 
curves such as those in Fig. 8 resulted. Note that 
la used as the abscissa In this figure 
because the basic theory underlying the gages shows 
the power loss per unit temperature rise to be lin- 
early proportional to this quantity. The range of 
the representative curves shown on Fig. 8 encom- 
passes all the gages and provides an Indication of 
the maximum variation between the gages. At the 
lower values of shear, conduction effects In the 
plastic substrate cause a departure from linearity 
(see the curve for wire 24), and this had to be 
accounted for In the calibration, especially since 
low v-jlues of shear were expected in the region of 
separated flow. 

Co. ■ ' irlson of Computed and Measured Flow Field 

i.n this section, comparisons will be made oc 
the omputed and measured flow-field characteristics 
for M “ 0.788, where shock-induced bound.* ry- layer 
sepiratlon occurs. In considering these comparisons. 
It mi St be recalled that the computations at this 
tins* ipply only to free-fllght conditions and may not 
do eiy agree with calculations employing Identical 
turbulence models but utilizing the contours of the 
upper and lower walls of the channel as boundary con- 
ditions. The comparisons, then, of the numerical 
solutions and the data must be viewed in a 
qualitative sense. 

In Fig. 9, the calculated and measured surface 
pressure distributions are compared. From experi- 
ence with Invlscld flow-field computations. It Is 
felt that the mismatch In the location of the shock 
wave shown In this figure is due In part to the dif- 
ferences between the boundary conditions of the cal- 
culatlonii ind those that actually existed in the 
channel with contoured walls designed for M ■ 0.773. 
01 course, the turbulence lodellng also has a larg< 
effect on the shock-wave location. The small Mach 
ar. < Reynolds number differences between the 



Fig. 9 Comparison of ev''erlmental and computed 
pressure distribution on airfoil with 
shock-induced separation 


calcuLitionn and the experimental conditions, how- 
ever, are believed to he Insignificant. This ludp- 
nent Is based on a comparison of a series of conr 
putrd results with a specific turbulence moriel and 
frae-fllght boundary conditions for a range of Mach 
and Reynolds numbers.^ A nualltatlve conclusion of 
the comparison between the data and computations In 
Fig, 9 Is that the turbulence model employed In the 
present calculations fsee Computational Method) 
resulta In too high a nressure recovery, rhvslcallv, 
this would imply that the predicted thickness of the 
separated region Is too narrow. An eddv dlffuslvltv 
that la too large overall could cause the nar*‘ow 
separation region,^ For turbulence In equilibrium 
with the mean flow, there la a contradiction here In 
that thinner regions of viscous flow usually result 
In smaller values of eddv dlffuslvltv in the outer 
regions since r » Thla Implies that the chief 

turbulence modeling difficulties He with the model- 
ing of the inner regions, regions I and III of the 
sketch In the section on Computational Method. An 
alternative meana of achieving lower effective eddv 
dif fuslvltles in the modeling Involves use of the 
concept of rurbulence lag (e.g.. Refs. 4 and 9). 

Figure 10 shown the distribution of the local 
shear stress coefficients, based on the free-stream 
dynamic pressure, that were obtained from the skin- 
friction gages at Re^ ■ 10,7 » lO*" and 14 3 » 10^, 
For values of x/c less than 0.675, that in, ahead 
of separation, predictions could he based on the 
Marvin-Sheaf fer boundarv-laver program and the 
experimental pressure distributions. The difference 
between these lines Is an Indication of the magni- 
tude of the effect Reynolds number has on the local 
skln-frlctlon coefficient for tne attached flow. 

The data obtained with the akln-frlctlon gages 
are shirm in Fig. 10 for Re^, ■ 14.3 » 10® and fo.- 
RCp “ 10.7 * 10®, It Is noted that many of the wire 
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Fig, 10 nistrlhi tlon of skln-f rlctlop coefficient 
In the X resence of bounda.r-laver 
separation 


gAReii Indicated on Fig. 2 and in Table 1 are nut rep- 
resented on Fig. in. The reason for this Is that 
after calibration several of the ga<’eH broke durlnr 
the running of tlie experl'oe.u . The gage vlre i .ir> 
extremely delicate and were danwiged by particles of 
dirt In the air that were accidentally Introduced 
Into the channel. Fortunataly, the damaged wires 
were limited to values of x/c less than n.S. 

Usually, the wires were damaged by being broken or 
pulled away from their lend-ln posts. A few wires 
were completely pulled off the model. Visual exam- 
ination of Che broken wires before and after runs 
Indicated that no fragments of wire extended Into 
the airflow to .'suse ma.ior disturbances. Fortunately, 
the objectives of the present experiment tfere Identi- 
fied with processes taking place aft of the 
x/c >0.5 station, where none of the wires was 
broken. The Increased surface roughness probablv 
Introduced by the upstream broken wires Is not 
expected Co have affected the results beyond x/c > 
0.67S because of the rather thick viscous layeis 
there and the relatively long distance between the 
last o' the broken wires and the region of gages that 
remain id Intact. 

ji comparison of the experimental data at 
x/c " 0.276, 0.308, and 0.531, where the flow Is 
attached, with the Marvin- Sheaf fer calculations 
(Identical with those used Co "calibrate" the wire 
gages) is an Indication of the repeatability of the 
skln-f rlctlon measurements. These points are all 
within lOX of Che boundary- layer computations corre- 
sponding to the appropriate Reynolds number and 
experimental pressure distribution. If the repeat- 
ability Is taken as a sKasure of the uncertainty 'f 
the measurement and If the percentage uncertainty In 
the power iseasurement is Created as constant, the 
uncertainty in skin friction at the low values In 
Che separated region Is estimated to be about 15Z. 

The uncertainty In the region of separated flow and 
low skin friction Is larger because the heat trans- 
ferred to the subotrace material of the heated wire 
gage becomes proportionately larger. 

The solid line represents the values of the 
skln-frlctlon coefficient obtained from Delwert's 
Navier-Scokes code and the turbulence modeling 
described previously for a Reynolds number of 
10 X 10^, and is Co be compared with the experimental 
data represented by circles. Again, because of the 
difference In boundary conditions between the free- 
fllghc cosgiutatlons and those that existed in the 
channel, only qualitative comparlsi as are meaningful. 
As with the pressure distributions, the predicted 
skln-frlctlon coefficient does not drop at quite the 
same chordwlae station as Che data. The predicted 
values indicate t'le shock-wave location at a larger 
x/c than the data, which Is generally consistent 
with Che pressure distribution results; however, the 
difference In shock location Indlcatef hy the skin- 
friction data Is sosmt/hat less than In the case of 
the pressure data. 

In general, Che experimental data and the cois- 
putations show a similar character in the separated 
flow region. After sepa-atlon, a maximum shear in 
the negative direction develops. This shear dimin- 
ishes in magnitude with distance back on the airfoil 
for a while, and Chen approaches a second maxlmvtm In 
the absolute sense as the trslling edge Is approached. 
This general character suggests strong flow reversal. 
Tlie computations tonflne this behavior to a region 
closer to the trslling edge Chun is evlden' in 'he 
da's, which ie conslscent with the downstream 


Predict ion of the shock-wave poaitlon. The mnanl- 
tiidc of the reversed shears, however, are aerlouslv 
imdernredicted Just he/rnd separ.ttion and near the 
iraMlnp edge. Near the trailing edge tliis miv be 
caused by a reversal in the pressure gradient in the 
computations bevond x/c ■ 0,9. It should be noted 
that bevond the airfoil the pressure gradVnt 
becomes strongly adverse again, which suggests a 
computational singularity at the trailing edge. 

Fluctuating Skl n -Prictlon Me. ’surements 

To gain some knowledge or the relative rharar- 
terlstlcs of skln-frlctlon fluctuations in regions 
of attached and separated flows, the ac coomonents 
of the outout voltages of the heated wires were 
recorded on magnetic tape for aubseouent analysis. 
Because the wire gages are imbedded in plastic, 
their frequency response Is low (1200 to 1600 Ha) 
compared to that of a conventional hot-wire anemom- 
eter; hence, the fluctuating voltage output hy the 
gages is In response to onlv some fraction of the 
total fluctuation In skin friction present at the 
airfoil surface. In order to compare the fluctuat- 
ing akin friction on an equitable basla of frequency 
content, the recorded signal of each heated wire was 
passed through a 1000-Hz cutoff filter prior to 
evalustlnf; Its rms voltage. The rms valuer of 
(liyPw'^w) were evaluated using the rms voltage output , 
AT, the resistance of each wire gage, and aasuning 
that the power loss per unit temperature varies 
linearly with (uyOyT„)*'^ over the range of the 
fluctuation. 

The chordwlse distribution over the airfoil of 
the ratio of the rms-to-mean vali'es of (tiv i" 

ahovn In ®'lg. 11 for shock- Ind-.ced separation at 
M ■ 0.7R4 and Re^ ■ P » 10". j.hcad ot the shock 
wave, in the region of attached flow, the fluctua- 
tions are a small fraction of the mean value and 



x/c 

"It. 11 Skln-frlctlon gage fluctuatlona (0-1000 II*) 
for thock-lndu'-ed separation on the airfoil; 

M - 0.784. Fe - 8 • 10‘ 

- c,- 



inrreane downMtream na the ahock la approached. In 
Che aeparaCed flow region behind the ahock wav^-<, 
where the nean value of la arwll ,'aer Tl'. l('), 

fluctuatlona in the quantity (i>w‘''w'w^ appro i h 
valiiea Chat are nearly aa large aa the mean and, near 
the tralliiig edge, again become a amall fraction 
thereof. The aign of Che ordinate In Fig. 11 rorre- 
aponda to that of In Che mean value and waa 

determined by the algn of tite akln-f rlctlon coeffi- 
cient In Fig. 10. It appeara Chat the fluctuatlona 
in "akin friction" at 1000 Hr or leaa are a large 
fraction jf the mean value In reglona of aeparated 
flow. 

Concluding Remarka 

The flush hot-wire akln-frlctlon gage has been 
demonstrated t-j be able to measure the local akin 
friction on a transonic airfoil experiencing shock- 
tnve-lndmed turbulenc-boundarv-layer separation In 
high Reynolds number flow. The estimated accuracy 
of these measarements la about ISt. Together with 
local surface pressures measured on the same airfoil, 
these data provide the bases for testing and verify- 
ing conputatlona based on RiynolHs-averaged Navier- 
SCukes equations that are closed through turbulence 
modeling. To date, however, the computer codes have 
not orogressed sufficiently to permit efficient 
development of a new turbulence model. One serious 
shortcoming of computer codes employed to date is 
th It they appl" strictly to frec-fllght conditions 
and do not accommodate test channel wall configura- 
tions. This mismatch In boundary conditions can 
have significant effects on the flow Held Immedi- 
ately next to the airfoil surface. Therefore, before 
different turbulence models can be quanLltatively 
tested against the pressure and skln-frlctlon dat.i It 
Is necessary to incorporate the channel wall config- 
uration and the boundary-layer displacement effects 
Into the boundary condition of tlie computer codes. 
This Is being done currently. 

With the proper boundary conditions. It will be 
possible to test computer codes containing different 
kinds of turbulence models and numerical logic 
against the data presented here. The predicted nres- 
sures already suggest that the turbulence model used 
in the present calculations yields a separation 
region that la narrower than would he consistent 
with the pressure measurements. The message wlcMn 
the akln-frlctlon data Is less cleat because the 
flow In the separation bubble Is senbitive to iiu'>del- 
Ing changes throughout the flow field end must await 
systematic turbulence model modif^c.ttlons In a code 
with proper boundary condltloas to provide modeling 
guidance. 

The comewncs above Imply tlie belief that 
Reynolds stress modeling can ultimately be made to 
fit both the experimental pressure and skln-frlctlon 
data In the separated layer. In the narrow enqilrl- 
cal sense this Is true, yet the large dynamic char v - 
ter of the skln-frlctlon data found to exist In the 
separated zone may be Incompatible with universal 
statistical turbulence models appropriate to other 
kinds of flow that are steady. 
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